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Abstract. Diffusion in a one dimensional random force field leads to interesting 
localisation effects, which we study using the equivalence with a directed walk model 
with traps. We show that although the average dispersion of positions (x 2 ) — (x) 2 
diverges for long times, the probability that two particles occupy the same site tends 
to a finite constant in the small bias phase of the model. Interestingly, the long time 
properties of this off-equilibrium, aging phase is similar to the equilibrium phase of the 
Random Energy Model. 



1. Introduction 

The properties of random walks in random environments can be markedly different from 
those of homogeneous random walks KJ. For example, the typical distance travelled by 
a diffusing particle in an unbiased random force field in one dimension grows with 
time as x oc log 2 t, instead of the usual \/t law 0. This is due to the fact that the 
potential energy typically grows as y/x, leading to very high barriers which slow down 
the progression of the particle. More strikingly, Golosov has shown that the relative 
distance between two particles in the same random force field remains finite even for 
large times , whereas it also grows as \/i in a homogeneous medium. This remarkable 
classical localisation phenomenon is due to the fact that the 'best' potential minimum 
which can be reached by the particles after a long time t is so much better than the 
'second best' that all the particles have time to gather there, before eventually moving 
to a still better location. 

In the presence of a non-zero average bias Fq > 0, several regimes must still be 
distinguished. For small enough F , the mean position of the particles grows as t^, 
where the exponent /x < 1 is proportional to F [[]], |, ^ []]. Beyond a critical force, the 
particles move to the right with a non-zero velocity. However, the dispersion around 
the mean velocity is still anomalous until /x reaches the value /x = 2. Beyond /x = 2, the 
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spreading is 'normal', i.e. Gaussian with a width growing as y/i. The question we want 
to address in this paper is whether the Golosov phenomenon survives in the presence of 
a non-zero average force. We will actually show that different 'localisation' criteria lead 
to different answers: while the average width of a packet of particles diverges with time, 
there is a finite probability (even at long times) that two particles are nearby in space. 
The physical picture is that the density of particles is concentrated on a finite number 
of sites, but the relative distance between these peaks grows with time. As we shall also 
discuss, there is a strong analogy between this problem and the low temperature phase 
of Derrida's Random Energy Model |7| . 



2. Model and simulation 



Actually, one can map the long time behaviour of the problem onto that of a much 
simpler directed walk model |1|, [], §, 0, where each particle hops to the right on an 
even-spaced discrete lattice, with a site-dependent hopping rate W n distributed as: 



For < 1 the average trapping time 1/W is infinite, this leads to the anomalous 
behaviour of the particles average position reported above. The probability to find the 
particle on site n obeys the following Master equation: 

dP (t) 

= -W n P n (t) + W n ^P n ^(t) (2) 

The properties of equation (H) can be discussed in detail || [K| . In particular, one can 
compute the average position of the particle, defined as: 



(x(t)) = J2nP n (t) (3) 



n=0 



where the overbar denotes the average over the W's. For // < 1, one easily finds 
(x(t)) = sm(iifi)/[7ifiT(fi + at large times 0. One can also compute the average 
width A 2 of the packet, defined as: 



A 2 = (x 2 ) - (x> 2 = £ n*P n (t) - £ nP n (t) (4) 

n=0 \n=0 / 

As shown in [[|, the width grows to infinity as 

A 2 = C{n)t 2 ^ (5) 

where C(/z) is a certain fi dependent number which can be explicitely computed ||, and 
which goes to zero for fi = 0. Naively, this means that the particles' relative positions 
get further and further apart (for // > 0) as t becomes large, at variance with Golosov's 



result for the unbiased case, which shows that (x 2 ) — (x) 2 remains finite for large t. 
However, one can still ask the following question: what is the total probability that two 
particles, initially at site n = 0, occupy the same site after time t ? This is obtained as 

oo 

Y 2 {t) = Y.Pnitf (6) 

n=0 

The notation Y 2 is introduced in analogy with spin-glasses, where the same question is 
asked about two copies (replicas) of the same system in equilibrium, and it measures 
the probability that these two copies occupy the same state [[□]]. Note that Y 2 is also 
often taken to be an indicator of localisation in quantum problems, where P n = \ip n \ 2 is 



the quantum probability of presence |12| . One can actually study generalized objects, 
such as: 

00 

Y q (t) = £ Pnitf (7) 
n=0 

which measure the probability that q particles occupy the same site. We have studied 
numerically Y 2 and Y 3 for the problem defined by equation (Q). In figure [I], we show the 
average value Y 2 {t) as a function of t for \i — 0.4 < 1. This quantity clearly tends to a 




Figure 1. Time evolution of the quantity Y%{t) in the simulation of a directed random 
model with hopping rates distributed according to a gamma distribution of index 
/i = 0.4. The simulation was carried out in a 20.000 point lattice with 1.000 particles 
and averaged over 250 disorder samples. Up to the last observation times more than 
99% of the particles remained within the lattice. This curve shows that Y^it) tends to 
a non-zero constant for asymptotic times. 



non-zero constant y 2 ([i), which we plot as a function of /1 in figure 0, together with our 
theoretical prediction (see below). 
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Figure 2. y 2 versus \i as given by numerical simulations (small circles), the analytic 
formula obtained in the text (solid line) and by the equilibrium result of the associated 
Random Energy Model, yi = 1 — fj,, (dotted line). The agreement with the theoretical 
prediction is good, except for values near fj, = 1: the observation time being finite, we 
expect to overestimate 1/2 (/■*) m ore and more the closer one gets to fi — 1 since the 
approach to the actual asymptotic value becomes very slow (logarithmic) for fj, = 1. 

We see that y 2 (fj, 0) = 1, as expected from Golosov's results, while we observe 
the tendency t/ 2 (/i — > 1) — >• 0, although our numerical data is biased in this limit due to 
the slow approach to the actual asymptotic value. In figure 0, we show the parametric 
plot of t/3 (/x) [i.e. the asymptotic value of Y 3 (t)} versus 2/2 (/•*)• Interestingly, the resulting 
curve is seen to be very close to 7/3 = £/ 2 (l + 2/2) /2 obtained within the so-called 'one 
step replica symmetry breaking' solution of equilibrium random systems. 

Before giving a more detailed physical interpretation of these results, we first turn 
to an analytic calculation of y 2 (/x)- 



3. Analytic derivation of 2/2 (m) 

In the trapping model that we study, the time that the particle sojourns in the z-th 
trap is given by t; = u{Ti where Tj is the characteristic trapping time of the z-th trap 
(Tj = Wf 1 ) and Ui is an exponentially distributed variable accounting for the individual 
thermal behaviour of the particle. It is then straightforward to write the following 
equation for the probability of a particle being at site n at time t, given a realization of 
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Figure 3. Representation of 1/3 versus 2/2 a s obtained from our simulation (small 
circles) and as given by the relation 2/3 = 2/2(1 + 2/2 )/2, characteristic of replica 
symmetry breaking (solid line). The agreement is remarkable and hints to some 
deep relation between the present dynamical problem and the equilibrium phase of 
disordered systems. 



the disorder {r^} (or, equivalently, {W^}): 

„ n I n \ / n— 1 \ / n \ 

Pnit) = J J! dinexp ( - E Ujj [t - E M * r *J # (E M * r * - *J ( 8 ) 

being the Heaviside function. Equation @ says that P(n, t) is the sum of the 
probabilities of all the possible thermal histories such that the particle has already 
done n — 1 jumps up to time t but not yet n. 

We are interested in the probability that two particles, starting together at t — 
remain at the same site after a sufficiently long time. Therefore, the quantity of interest 
is the probability that two particles are at site n at time t: 



/ft, I It, \ I I V J- \ I Hi 

J] duidvi exp - E(% + vj) j f t - E J # I E 



n-1 



x^ft-E^^E^-*) ( 9 ) 

or, in the Laplace domain, 

P 2 n {E) = / II ditto* exp - + «,■) / ' 8 ' . e^d* (10) 
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Integrating and using 6 functions we can rewrite this last expression as 



„ n 

EP 2 n {E)=2 j \\dudvi 



n-l 



exp —E UiTi — exp I —E ^ 



/n-l 



n-l 



X 



+ 



n-l 



exp -E ^ UiTi - exp -E ^ Vin 



/n-l 



n-l 



X 



J2( U i ~ V i) T i ) \ J2 V i T i-J2 U i T i f eX P ~ J2( U 3 + 



where the three #'s in each summand implement the maximum condition in the lower 
limit of the integral, the minimum condition in the upper limit and the condition that 
the upper limit is greater than the lower limit in fllCf), respectively. In (11) we have also 
used the symmetrical integration with respect to Ui and i>j to simplify somewhat the 
expression (whence the factor 2). 

It can now be proved that expression (11) is equivalent to a much simpler 
formula, where only one 6 function appears per summand and this is accomplished by 
conveniently renaming Ui <-» V{ as integration variables for some terms and by bearing 
in mind that Ui, Vi, and are all positive, whence the summations are all monotonous 
increasing functions of the step number n. This procedure leads to 



„ n / n 

EP^(E) = 2 J J] duidvi exp I - £>i + 



n-l 



/n-l 



exp —E J2 U i T i )° \J2( U i ~ V i) T i 



n-l 



exp —E J2 U i T i J2( V i ~ U i) T i ) ~ eX P \~ E J2 U i T i 



/n-l 



/n-l 



UiTi 



(12) 



In order to proceed we now use the following representation of the 9 function in the 
complex plane 

e i>~ 



A 



-dX 



and (12) turns into 

EPl{E) = — ^ / J] d Ul dv % exp - 5>i + 

71 J— co A J \ 



(13) 



n — 1 



exp y~Y^[Eui - i\(ui - Vi)]TiJ 

n-l \ 

- exp ( - Y} Eu i - iX (vi - Ui))Ti ) - exp f - £ - iA(t*i - ^)h ) e - iAt, " r « 

(£— iA)u„r n 



n-l 



+ exp - X \ Eu i - l K v i - 



Ui m e 



(14) 



This formulation of the equation permits the factorization, within each summand, of 
the contributions of each duidvi so as to get 

raa J\ rn-1 

EP 2 n {E) 



i r°° dX 

7T J-oo A 



L> -L th lb n 

[] / dudve- u - v e- [Eu - iX{u - v)]n - JJ / dudve~ u ~ v e- [Eu ' iX{ - v ~ u 

i i 

/n— 1 „ 

i 

/n—1 . 
du n e- (1+ETn+iXTn)Un J| / dudve- u - v e- lEu -' lX{v ' u)]n 



If we now define the following functions 
F(E,X) = 



dudve~ u - 1 ' e -[Eu-i\(u-v)]n 
G(E, X) = J due~ u e-( Eu + iXu ^ 



(15) 

(16) 
(17) 



where the bar over the exponential stands for the average over the possible values of Tj, 
the disorder average of (15) is readily written as 



EP*(E) 



d\ 

. — {[1 - G(0, A)] F(E, X) n - [F(E, -A) - G(E, A)] F(E, -A)"} (18) 

7T J-oo A 



We now sum ([18]) for all the values of n in order to obtain the function Y^E): 



— i r°° dX 
Y 2 (E) = -f -y 

X3 A 



7T 



1-G(0,A) F(E,-\)-G(E,\) 



(19) 



_1-F(£,A) 1-F(£,-A) 

Since = W^ -1 the functions F(E, A) and G(E, A) can be expressed in terms of the 
distribution of hopping rates p{W): 

W W 



F(E,X) 



G(E,X) 



dWp(W) 



W + E - iXW + iX 



1 - 



A 2 



E - 2iA Jo 



dW 



p(W) (E-iXf 



W + iX E - 2iA Jo 



dW- 



p{W) 



W + E - iX 



dWp{W) 



W 



W + E + iX 



1 - (E + iX) [ 
Jo 



"m-. p(w) 



W + E + iX 

Using these expressions in ( |19"D we get, after some algebra, 



(20) 



(21) 



9 r°° 

EY 2 (E) = - I Re 



7T JO 



f(iEu)-f(E-iEu) 
u 2 f{iEu) + (1 - iu) 2 f(E - iEu) 



[l-iu] 



du (22) 



the function f(z) being given by the integral 

1 



/(*) 



dWp(W) 



W + z 
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The result ( |22| ) is now straightforwardly applied to the relevant distribution of 
hopping rates p{W). For instance, it is reassuring to see that for a non-disordered 
lattice, p(W) = 5{W - W ), equation (|22]) yields: 

Y 2 (t) ~ / 

which can be obtained directly since the two particles are independent, and this means 
that the probability that the two particles are on the same site is inversely proportional 
to the typical distance. 

For the case that we are exploring here, we focus on a distribution of the kind ([]]), 
whence the function f{Ez) turns out to be 



f{Ez) 



dW 



W n-x e -w 



1 



OO <£, 

dx 



o T(/j)x + z 



T(fi) W + Ez 

assuming in the last equality that E is sufficiently small so as to neglect the factor e~ Ex 
within the integral. Some further developments making use of the definition of gamma 
functions allow one to get to 



f{Ez) ~ T(l - fijE^z^ 1 as E 
Introducing this expression in (E22|) we finally obtain 



— , . 2 



EY 2 (E) 



Re 



(iuy- 1 - (1 - ru) M_1 







1 



du 



(23) 



7T Jo l(l — iu)^ +1 — (iu 

which is indeed a finite integral when < \i < 1. Equation (|23"D can now be trivially 
transformed to the time domain again and we obtain a //-dependent constant asymptotic 
result: 



vM = - 



IT JO 



Re 



111) 



m 



(111 



111 



du 



(24) 



We have calculated (^4|) numerically for different values of /x in the interval of interest 
and we have compared the results to the data obtained from the simulation in figure 
H]. The agreement is quite good except for values near \i — 1, where the results of the 
simulation turn less reliable because of the slow relaxation to the actual asymtotic value 
in (EE 



4. Discussion 



How can one reconcile the fact that, at the same time, the typical distance between 
two particles grows with time (as f 1 ) and that the probability to find them at the 
same site tends to a finite constant? The physical picture is that of figure £|, where 
the probability P n (t) is shown for a single sample and a fixed t. One sees that this 
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Figure 4. Distribution of probability after a time t = 7 x lO 10 !^ for a particular 
sample of disorder in our ID directed random walk model with (i = 0.4. The simulation 
was done with 1.000 particles in a lattice of 20.000 sites. 



probability distribution is made of several sharp peaks that gather a finite fraction of 
the particles. However, the position of these peaks is scattered on a region of space of 
width £ M . As time progresses, the position and relative weights of these peaks of course 
change, but at any given (large) time only a finite number of peaks, corresponding to 
very large trapping times, contain most of the particles. This is clearly related to the fact 
that, for a given particle, most of its life was spent in the deepest 'trap' encountered 



up to time t Hl3| , p~4| . This behaviour is typical of the Levy sums: when fi < 1 the 
sum of individual trapping times is dominated by the largest one. There is also a 
strong connection with the physics of the Random Energy Model in its low temperature 
(glassy) phase. The distribution of Boltzmann weights (which are also the residence 
times within each state) is there again a power-law with an exponent fi < 1, whence 
only a finite (but random) number of states contribute to the full partition function 
fl~5| , p!6| , and the probability that two independent copies of the same system occupy 
the same state is finite (and equal, on average, to yi = 1 — /i). As discussed in detail in 
fl"6|| , this is in turn related to 'replica symmetry breaking'. All the y q J s can be computed 
and one finds, in particular, y 3 = 7/2(2/2 + l)/2. As explained above, we find analytically 
that y 2 7^ 1 — /i in the dynamical model, which means that the system can never be 



considered in equilibrium, although the dynamics gets slower and slower with time [T7 
At the same time, however, the equilibrium relation 7/3 = 7/2(2/2 + l)/2 appears to be 
fulfilled (see figure 3), suggesting that some kind of pseudo-equilibrium can be defined, 
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for which equilibrium methods such as the replica method could be applied. It would 
be interesting to extend the method of the previous paragraph to calculate all the ?/ g 's 
exactly, and to check whether they agree with the replica prediction. 

Finally, it is interesting to note that the above biased model exhibits aging effects 
when fj, < 1 [[IB], [H| [TI| pPf . In this context, a classification of different aging models was 



proposed in [21], in terms of the asymptotic 'clone overlap' function. The idea is to look 
at the evolution of two identical systems (replicas), driven by the same thermal noise 
until t = t w , and by independent thermal noise for later times. The two replicas can 
either separate with time (type I aging), or remain close even after infinite time (type 
II aging). One sees from the above example that, depending on the way in which one 
measures the 'closeness' of the two particles, one concludes differently. This situation 
is reminiscent of the quantum localisation model introduced in |[22|| , where states are 
(in certain regions of parameter space) both extended and localised, depending on the 
property which is studied. The physical nature of these quantum mixed states is actually 
very similar to the one discussed above. 
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